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Abstract 

Results for the isovector axial form factors of the proton from a lattice QCD 
calculation are presented for both point-split and local currents. They are ob- 
tained on a quenched 16 3 x 24 lattice at = 6.0 with Wilson fermions for a 
range of quark masses from strange to charm. For each quark mass, we find 
that the axial form factor falls off slower than the corresponding proton electric 
form factor. Results extrapolated to the chiral limit show that the q 2 depen- 
dence of the axial form factor agrees reasonably well with experiment. The 
axial coupling constant qa calculated for the local and the point-split currents 
is about 6% and 12% smaller than the experimental value respectively. 
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As have electromagnetic form factors, the isovector axial form factor qa{q 2 ) has 
been commonly used to constrain the construction of models of the nucleon in order 
to incorporate PCAC and the Goldberger-Treiman relation [fll 0]. Similar to vector 
meson dominance in electromagnetic form factors, the isovector axial form factor 
seems to be quite sensitive to whether degrees of freedom in the A\ channel (e.g. 
pir,uj7nr or Ai itself) are introduced in the effective theory || [|, [|. In view of the 
fact that the EM form factors and the magnetic moments of the nucleon in recent 
lattice QCD calculations are within 10 to 15% of the experimental results [0, ^J, it is 
natural to extend the study to the axial form factor and to determine to what extent 
the experimental results can be reproduced in lattice QCD, subject to the limitations 
of the quenched approximation, finite size effects and the extrapolation to the chiral 
limit. Especially interesting is to check the q 2 dependence of the axial vs. the electric 
form factors to see if the experimental difference between them is borne out in the 
lattice calculation; the q 2 dependence of their ratio should be fairly independent of 
the systematic errors due to the present limitations of the lattice calculation. In 
this paper, we study the axial form factor in lattice QCD with Wilson fermions and 
compare it to the electric form factor calculated previously and to experiment. We 
examine its systematics as a function of the quark mass from strange to twice that of 
the charm mass, and compare the axial coupling constant qa to previous quenched 
calculations with different volumes || ||. [il], |llj. Fixing qa to the non-relativistic value 
of 5/3, we determine the finite lattice renormalization for heavy quarks. Assuming 
axial dominance, we extract the AiNN form factor. 

Lattice gauge calculations have been carried out to study the electromagnetic form 
factors of the pion[|l^] and the nucleon|| [7|]. The same sequential source technique 



(SST) using the zero-momentum point nucleon interpolating field as the secondary 
source is applied here to study the axial form factors ]7) . In choosing the final hadron 
as the secondary source, one can sew the quark propagators together at the point 
where the current couples in the three-point function. This has the advantage of 
being able to study different currents at various momentum transfers. 

The lattice two- and three-point functions that we calculate are the following: 

G a P a p(t,P) =T / e^ s (0\T(x a (x)x Q (0)\0), (1) 

X 

Gf AP (tf,P,t,q) = E e- i ^ + ^(0|T(x a (^)^(^)x /3 (0))iO), (2) 

x f ,x 

where \ a is the proton interpolating field and A^(x) is either the point-split axial 
vector current 

A^ s = i2K[^(x)- lfll5 U^x)^(x + fi)+$(x + fi)-^ 5 Ul(x)^(x)}, (3) 



or the local current 



A L ° C = i2Knl>(x)wtfl>(x). (4) 
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Phenomenologically, the axial vector current matrix element is written as 



(ps|A*(°)Ip" s ') = u & s )(WA(g 2 ) - q^h A (q 2 )]^u(p' , s'). (5) 



It has been shown |K| that when tj — t and t » a, the lattice spacing, the ratio of 
eqs. ([|) and (|l|) gives the following relations for the A 3 and A 4 matrix elements 

G P a P (t f ,Q) 2E q [yMq ' E p + m Mq >h KJ 

T^Gf A4P (t f ,6,t,q) q 3 



e-^- m »ig A (q 2 ) + (m - E q )h A (q% (7) 



GF P {t f ,0) 2E, 



where T = ^ ^ 3 jj J . When <f = 0, eq. (|6|) reduces to the coupling con- 
stant. For finite momentum transfer, g A (q 2 ) can be obtained from eq. (|6]) upon 
eliminating the kinematic factor extracted from the ratio of the two-point func- 
tions Gpp(t,q)/Gp p (t,0), and by setting q 3 to zero. The only exception is when 
q 2 = 3(27r/La) 2 where q 3 can not be set to zero. In this case, we use both eqs. (|) 
and (0) to obtain g A (q 2 ) at this q 2 . 

Our quenched gauge configurations were generated on a 16 3 x 24 lattice at f3 = 6.0. 
The gauge field was thermalized for 5000 sweeps from a cold start and 12 configu- 
rations separated by at least 1000 sweeps were used. Periodic boundary conditions 
were imposed on the quark fields in the spatial directions. In the time direction, the 
fixed (uncoupled) quark boundary condition was used. All quark propagators were 
chosen to originate from lattice time slice 5; the secondary nucleon source was fixed 
at time slice 20 (except for k = 0.154 where the quark propagators from time slice 3 
to 22 are used). We also averaged over the directions of equivalent lattice momenta 
in each configuration; this reduces error bars. 

The results presented here are for Wilson fermions with k = 0.154, 0.152, 0.148, 
0.140, 0.133, 0.120 and 0.105 so that a range of quark masses from the strange to 
about twice the charm mass is covered. To extract the lattice axial charge g^, in- 
stead of using the ratio in eq. (|6]), we fit the three-point and two-point functions to 
two exponentials in the form of g A fe~ mt and fe~ mt simultaneously using the data- 
covariance matrix to account for the fact that they are measured on the same set of 



gauge configurations ||14j| . The range of t for the two-point functions are chosen to 
overlap with the tf in the three-point function. This is done except for the heavy 
masses (i.e. k = 0.120 and 0.105 for the point-split current and 0.140 and 0.133 in 
addition for the local current) where eq. (0) is used. For q 2 ^ 0, we used the combined 
ratios in eqs. (|) and (0) and G PP (t,q)/G P c p(t,0) to extract g A (q 2 ), since we do not 
have enough gauge configurations to warrant a simultaneous fit ||15|| . The errors are 
obtained through the jackknife in this case. 
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Table 1: Unrenormalized coupling constants g\ calculated with the point-split and 
the local currents for different quark masses. Also listed are the \ 2 P er degree of 
freedom from the covariance fit. 





.154 


.152 


.148 


.140 


.133 


.120 


.105 


(point-split) 


1.28(17) 


1.26(6) 


1.24(4) 


1.11(1) 


0.98(1) 


0.763(7) 


0.546(7) 


X 2 /N DF 


0.28 


0.03 


0.77 


0.43 


1.3 






g\ (local) 


1.47(18) 


1.43(5) 


1.40(5) 


1.25(1) 


1.11(1) 


0.88(1) 


0.651(6) 


X 2 /N DF 


0.20 


0.70 


0.77 











In Table 1, we list the unrenormalized coupling constants g\ = g^(q 2 = 0) for 
the point-split current (P-S.C.) and the local current (L.C.) for different k. We see 
that the results from the P-S.C. are lower than those of the L.C. We will address this 
point later. 

Plotted in Fig. 1 are the unrenormalized lattice isovector axial form factors g^iq 2 ) 
of the proton as a function of the quark mass in dimensionless units m q a = ln[l + 
1/2(1/k — l/ft c )] for different momentum transfers q 2 from to 4 times [2k /La) 2 for 
the point-split current. (N.B. q 2 = (E — m^) 2 — q 2 for the four-momentum transfer 
squared.) Also included are earlier results for g\ from smaller lattices (10 3 x 20 [0 
and 12 3 x 22). Although error bars overlap, we note that the finite volume effect is 
still appreciable for the light quarks. A recent quenched calculation with the same 
spatial dimension (16 3 x 40 with /3 = 6.0) |l tj shows that g\ at large k values of 0.154 
and 0.155 for the local current are at 1.36(4) and 1.39(14). Our result of 1.47(18) 
for k = 0.154 is compatible with these and is also in agreement with the earlier 
calculation with a renormalization improved Wilson action |TD[ . The extrapolation of 



g\ to the chiral limit at k c is carried out with the correlated fit to a linear dependence 
on the quark mass m q a for k = 0.154, 0.152 and 0.148. The covariance matrix is 
calculated with the single elimination jackknife error for g\, itself calculated from the 
simultaneous fit described above for different k. 

This extrapolation gives the unrenormalized g\ at the chiral limit to be 1.28(9) 
for the P-S.C. and 1.47(9) for the L.C. with x 2 / n df = 0.01 for both cases. Using the 



finite lattice renormalization Zi to be 0.86 for the P-S.C. [I6"| and 0.8 for the L.C. 
fPTfl as determined from current algebra matrix elements, the continuum which 
equals Z L g\ is 1.10(8) for the P-S.C. and 1.18(7) for the local current. This is about 
12% and 6% smaller than the experimental value of 1.254(6) respectively. 

It was demonstrated in a previous study |1 that the lattice axial charge g\ for 
the Wilson action is (5/3)2/t for the L.C. and ~ k 2 for the P-S.C. at the static limit. 
A quark mass dependent factor f(m q a) = e mqa = 1 + 1/2(1/ k — l/n c ) has been 
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Table 2: Renormalization factor Z\ for the heavy quarks 



K 


.140 


.133 


.120 


.105 


Z L A (local) 


0.98(1) 


0.97(1) 


0.97(1) 


1.00(1) 


Z A (point-split) 


1.10(1) 


1.09(1) 


1.11(1) 


1.19(2) 



introduced to local scalar [TT| and axial currents flHj, [TT| to smoothly interpolate 
between the limits of m q — > where / — > 1 and m q — > oo where / — > I /(2k) so 
that the static limit from the Wilson action matches with the continuum results. A 
similar factor is found in analyzing the Wilson action for non-relativistic quarks []18[ . 
We multiply g\ by this factor e mqa for both the L.C and P-S.C. and plot them in 
Fig. 2 as a function of the quark mass m q a. It is seen that from k = 0.140 on, 
the curve flattens and the value is very close to the non-relativistic quark model 
result of 5/3. We believe that this means the onset of the non-relativistic limit 
is near k = 0.140 and additional fine tuning beyond the mean field finite lattice 
renormalization factor e mqCL is needed to bring g& to be 5/3 for quark masses heavier 
than k = 0.140. In other words, the lattice field if>L is related to the continuum field 
ip c through ip c = (2Ke m <> a Z%) l / 2 i) L . The factor Z\ is listed in Table 2 for both the 
L.C. and the P-S.C. for the heavy quarks. To stress the fact that these k's are not at 
the static limit, we should mention that (local) for k = 0.105 is 2k(3.10) not the 
static value of 2k(5/3). It is the Z\e mqa factor that brings it to the non-relativistic 
value of 5/3. It is worthwhile pointing out that the finite lattice renormalization can 
be done for heavy quarks in a gauge invariant way without having to renormalize to 
quark properties which involves gauge fixing. Z\ for these quark masses is found to 
be very close to unity for the L.C. This means that the 2Ke mqa is a satisfactory ansatz 
for the renormalization. Z\ for the P-S.C. turns out to be 10 to 20% larger. One 
would think that this is due to the extra gauge-field link in the P-S.C. (eq. (HP) and 
that replacing it with the mean field average Uq which is between 0.8 and 0.9 at (3 = 6 
would reconcile the difference between the Z\ 's in these two currents. However, as 
we pointed out earlier, the g\ for the P-S.C. falls faster than that for the L.C. by an 
extra n toward the static limit ||. Hence, Z\ = Z A (n,g) is a function of both k and 
g. It can not be related to the mean field alone. 

For q 2 ^ 0, we extrapolate the form factor to the chiral limit with the procedure 
described in Ref. [|TTJ. The error of the linear extrapolation is calculated from the 
single elimination jackknife method. For each jackknife sample, the fit is carried out 
with the covariance matrix to calculate y 2 /N DF . The average x 2, /N DF turns out to 
be about 1 or less. 

In Fig. 3, we plot the axial form factor extrapolated to the chiral limit for the L.C. 
and the P-S.C. in comparison with the experimental result. In doing so, we have used 
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the calculated nucleon mass to set the scale for the momentum as was done in ref . [[7] . 
For comparison, we also show the calculated electric form factor Ge extrapolated to 
the chiral limit [|7] (without the covariance matrix) and the corresponding experimen- 
tal results. The experimental (m(<? 2 ) has been measured in neutrino-neutron scattering 
and pion electroproduction. The neutrino data gives a good fit in the dipole form 
up to q 2 = 3GeV 2 /c 2 @, i.e. g A (q 2 ) = g A (0)/(l - q 2 /M 2 A ) 2 , with the axial vector 
coupling constant g A (0) = 1.254±0.006 and M A = 1.032±0.036GeV (world average). 
Our fit of the axial form factor to the dipole form yields M A = 1.03±0.05GeV for the 
P-S.C. and M A = 1.03±0.03GeV for the L.C. which are very close to the experimental 
dipole mass. Similarly, the fitted dipole mass for Ge is Me = 0.77±0.03GeV which is 
close to the experimental dipole mass of 0.828 ± 0.006GeV. We should quickly point 
out that the favorable agreement with the experiments to a few per cent level for the 
dipole mass could be fortuitous at this stage since we have not adequately assessed 
the systematic errors due to the finite size effect, the dynamical fermion effect, etc. 
Especially, the isovector part of the charge radius of the nucleon has a chiral him 2 
divergence while approaching the chiral limit P0|. This correction to the lattice ex- 



trapolation to the physical pion mass can increase the charge radius by a significant 
amount However, unlike the isovector vector current which can couple to two 
7r's leading to a pion loop, the isovector axial current couples to three 7r's. Hence, it 
is not subjected to the chiral Inm 2 correction as alluded to for the charge form fac- 
tor. Notwithstanding these corrections (e.g. finite volume, dynamical fermion, chiral 
limit), we want to stress the fact that the q 2 dependence of the axial form factor 
shows that it is harder (i.e. falls off slower) than the corresponding GE{q 2 ) at each 
quark mass we studied. This feature, shown at the chiral limit in Fig. 3, is consistent 
with the experimental data (i.e. the dipole mass of 1.03 GeV for the axial case is 
higher than the 0.828 GeV for Ge)- This feature, we believe, is likely to survive the 
various lattice and chiral corrections. 

As far as the q 2 dependence is concerned, we find that the falloff is faster as the 
quark mass decreases. This is evidenced in Fig. 4(a) where g A {q 2 )/g A {0) for the 
P-S.C. is plotted for different quark masses. This is consistent with the fact that the 
associated meson cloud, specifically the A\, will have a larger Compton wavelength as 
the meson masses decrease with the quark mass. In view of the success of the vector 
dominance and the Goldberger-Treiman relation, the idea of axial dominance for the 
isovector axial form factor seems to work well phenomenologically []5J . With the axial 
dominance assumption in mind, gA{q 2 ) can be written as 

g A (q 2 ) ~ g Al NN(q 2 )/(l - q 2 /m 2 M ) (8) 

where g Al NN(<I 2 ) is the A 1 NN form factor. This is calculated by g A (q 2 )(l — q 2 /m Ai ) 
at k = 0.148, 0.152, and 0.154 and extrapolated to the chiral limit. The results for the 
P-S.C. normalized at q 2 = 0, are plotted in Fig. 4(b). We see that the q 2 falloff is now 
less sensitive to the quark mass than the case of g A (q 2 )- Finally we fit the gA 1 NN(q 2 ) 
obtained this way with a monopole and found a monopole mass of 1.08 ± 0.06GeV 



5 



for the P-S.C. and 1.01 ± 0.04GeV for the L.C. 

To conclude, we have calculated the isovector axial form factor of the nucleon for 
quark masses from strange to two times the charm mass. Albeit it is a quenched 
approximation, we find that g^, extrapolated to the chiral limit, is about 6% smaller 
than in experiment for the local current and about 12% smaller for the point-split 
current. It is important to include the e mqa correction factor for heavy quarks. The 
additional finite lattice correction for the local and point-split axial currents has been 
determined for heavy quarks in a gauge invariant way. We note that <m(<? 2 ) is harder 
than Ge{(1 2 ) of the proton which is in agreement with experiment. This finding is 
most likely to be preserved when finite volume, chiral, and dynamical fermion effects 
are included. Assuming axial vector dominance, we have extracted the A\NN form 
factor. For future studies, it is essential to improve the calculation by expanding the 
volume and incorporating dynamical fermions. 

This work is partially supported by the DOE Grand Challenge Award, DOE Grant 
No. DE-FG05-84ER40154 and NSF Grants No. STI-9108764 and PHY- 9203306. We 
would like to thank R.M. Woloshyn for providing us with his unpublished results on 
the 12 3 x 22 lattice. 
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Figure Captions 

Fig. 1 The unrenormalized isovector axial form factor g\(q 2 ) for the point-split 
current as a function of the quark mass m q a. The top curve is for q 2 = 0, the rest are 
for q 2 from 1 to 4 times of (2tt/ La) 2 in descending order. 

Fig. 2 The lattice g\ multiplied by the finite lattice correction factor e mqa for the 
point-split and the local currents as a function of the quark mass. 

Fig. 3 The calculated 5u(<? 2 ) for the L.C. and Gsiq 2 ) of the proton as a function 
of — q 2 . The dashed line is the fit to the experimental gA{q 2 ) with a dipole mass of 
1.03GeV. The solid curve is the same for the experimental GE{q 2 ) with a dipole mass 
of 0.828GeV. 

Fig. 4 (a) g^iq 2 )/ '9a f° r the P-S.C. for different quark masses as a function of —q 2 a 2 . 
(b) gAxNN^q 2 ) I gA 1 NN{0) defined in eq. (JD from the P-S.C. for several quark masses. 
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